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ABSTRACT 

This  study  evaluates  the  accuracy  of  an  established 
reliability  measurement  procedure  (NAVWEPS  OD  29304)  by 
computer  simulation.   The  reliability  measurement  procedure 
assumes  components  fail  according  to  an  Exponential  Failure 
Law.   This  study  tests  the  accuracy  of  that  procedure  when 
components  obey  a  rfelbull  Failure  Law  or  a  Log  Normal  Fail- 
ure Law. 
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CHAPTER  I 
INTRODUCTION 

The  highly  complex  weapon  systems  being  developed 
for  the  Navy  today  require  a  rapid  assessment  of  system  and 
subsystem  reliability  during  research  and  development,  and 
production  phases.   Currently  there  exists  a  statistical 
model,  the  Guide  Manual  For  Reliability  Measurement  program 
(NAVrtEPS  OD  29304),  that  "...can  be  utilized  by  all  contract- 
ors  for  subsystem  reliability  measurement  and  by  the  Navy 
for  weapon  system  reliability  measurement.  "1   This  model  is 
versatile  in  that  it  permits  the  combination  of  test  data 
from  all  levels  and  can  be  continuously  updated  as  new  test 
data  becomes  available .   The  model  is  thus  a  rapid  approx- 
imation procedure  for  determining  system  reliability.   How- 
ever, the  procedure  has  some  restrictive  assumptions  which 
must  be  kept  in  mind. 

One  notable  restrictive  assumption  is  that  all  compo- 
nents are  assumed  to  have  a  constant  failure  rate.   This  may 
not  be  necessarily  true,  although  many  times  is  a  good 
approximation,,   This  leads  to  the  commonly  used  Exponential 
Failure  Law0   However,  suppose  the  failures  were  actually 
given  by  a  ^eibull  or  a  Log  Normal  Failure  Law.   Then  this 
procedure  must  adequately  estimate  the  system  or  subsystem 


•'•Guide  Manual  For  Reliability  Measurement  program 
(NAVWEPS  OD  29304,  15  May  1965),  p.  3-1. 


reliability  under  the  Exponential  Law  assumption. 

The  purpose  of  this  study  is  to  evaluate  the  accuracy 
of  the  procedure  in  NAVtfEPS  OD  2930^  when  failures  actually 
occur  in  a  ^Jeibull  and  Log  Normal  fashion. 

The  approximation  procedure  uses  statistical  estimates 
of  failure  rates  based  on  sample  data  and  are  thus  subject 
to  statistical  uncertainty.   Therefore,  at  best,  the  pro-  f 
cedure  must  yield  a  lower  confidence  limit  on  system  or 
subsystem  reliability,  or  an  upper  limit  on  failure  rate. 
It  is  desired  to  test  the  accuracy  of  this  procedure  by 
simulating  the  distribution  of  the  lower  confidence  limit 
for  various  systems  and  comparing  the  simulation  results 
with  the  true  reliability. 

The  method  of  simulation  involves  obtaining  random 

Weibull  and  Log  Normal  failure  time  variates  for  components 

in  a  given  system  with  a  true  system  reliability,  R  . 

s 

These  quantities  represent  the  times  to  failure  of  compo- 
nents provided  they  are  less  than  the  duration  of  the  test 
(i.e.  the  planned  test  time).   Then,  using  the  statistical 
procedure  in  the  Guide  Manual  For  Reliability  Measurement 
Program,  a  100(1-  a )%   lower  confidence  limit,  R    ,  »,  for 
R  is  obtained.   A  large  number  of  replications  of  this 

process  are  completed  to  construct  a  distribution  of  R  T /  \ 

s ,L(a ) 

If  the  procedure  is  indeed  accurate,  the  (l-a)tjl  percentile 
point  of  the  constructed  distribution  will  be  R_.   That  is, 

S 

A 

if  Rs,L(ct)  ls  a  true  100(1 -a)^  lower  confidence  limit  for 
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*  ith 


Rs,    then  P(Rs>Rg   L(a\)   ■   1-a.      If   the    (l-a)™  percentile 

A 

point  of  the  distribution  of  Rg  L/a\  is  denoted  by  A-,    , 
then  An  _ a  should  equal  R  ,  regardless  of  the  set  of  param- 
eter values  chosen.   Thus  a  measure  of  the  accuracy  of  the 
procedure  is  given  by  the  quantity   A-,    -  fi  I  ,  Two  other 
measures  of  the  accuracy  of  this  procedure  are  the  mean  and 

A 

standard  deviation  of  the  distribution  of  Rs  T(a). 

The  accuracy  of  this  procedure  is  examined  for  a 
variety  of  combinations  of  components  and  sets  of  parameter 
values.   In  all  cases  the  systems  under  consideration  con- 
sist of  four  components  in  series.   All  components  are 
assumed  to  behave  independently  so  the  system  reliability 
is  the  product  of  the  four  component  reliabilities. 

As  expected,  the  closer  the  itfeibull  and  Log  Normal 
distributions  approximate  the  Exponential,  the  better  the 
accuracy  of  the  procedure „   However,  for  some  components 
with  failure  time  distributions  differing  considerably  from 
the  Exponential,  the  accuracy  is  fairly  good  using  this 
method  of  reliability  measurement.   If  the  absolute  differ- 


is  less  than  .02,  the  accuracy  is  called 
good.   The  accuracy  of  the  procedure  is  good  for  planned 
test  times  of  short  duration  and  diminishes  as  they  are 
increased.   The  study  concludes  that  the  reliability  measure- 
ment procedure,  being  an  approximation  procedure,  is  useful 
for  some  systems  following  tfeibull  and  Log  Normal  Failure 
Laws  provided  they  do  not  deviate  from  the  Exponential  Fail- 
ure Law  by  too  significant  an  amount. 


CHAPTER  II 


STATISTICAL  RELIABILITY  MODEL2 


The  statistical  Reliability  Model  provides  the  guide 
lines  for  reliability  measurement  on  Navy  weapons  systems. 
The  "Guide  Manual"  can  be  used  by  both  contractors  for  sub- 
system reliability  measurement  and  by  the  Navy  for  weapon 
system  reliability  measurement.   A  brief  description  of  this 
model's  assumptions  and  methodology  are  presented  in  this 
chapter. 

The  system  reliability  obtained  is  based  upon  testing 
components  for  a  time  duration,  called  the  planned  test  time, 
or  until  failure  occurs  under  a  given  stress  condition. 
Also,  the  model  provides  for  testing  a  system  for  component 
failure  by  the  same  method. 

This  model  considers  components  in  a  complex  system 
which  must  operate  successfully  over  a  defined  "mission", 
lasting  a  specified  duration  of  time.   Time  can  thus  be 
measured  in  "mission  units".   System  reliability  is  thus  a 
function  of  environmental  and  usage  stresses  which  include 
vibration,  shock,  etc.,  and  the  operating  vs.  non- ope rating 
conditions  as  well  as  a  function  of  the  associated  part 
failure  rate  parameters  together  with  the  time  duration  of 
the  environments  and  usage  stresses.   The  reliability 


2Gulde  Manual  For  Reliability  Measurement  Program 
(NAVWEPS  OD  29304,  15  May  1965),  pp.  3-1  to  3-15. 
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measurement  system  is  dependent  upon  the  following  assump- 
tions : 

1.  Constant  failure  rate.   The  exponential  failure 

law  is  assumed  to  hold,  mainly  for  its  common 
usage,  mathematical  simplicity,  practical 
simplicity,  and  reasonability. 

2.  Addltivity  of  stress  effects.   "The  failure  rate 

induced  by  two  simultaneously  acting  stresses 
is  equal  to  the  sum  of  the  failure  rates  due  to 
the  two  stresses  acting  sequentially. "3  Thus, 
mission  experience  can  be  simulated  by  adding 
data  from  separate  environmental  tests. 

3.  Independence  of  component  failures.   "...This  is 

assumed  because  components  are  normally  tested 
individually  by  type,  and  subsystem  reliability 
is  estimated  using  component  and  other  appli- 
cable test  results."^ 

4.  Failure  rate  constancy.   "The  failure  rate  is  con- 

sidered a  function  of  only  the  stress  acting. "5 
In  other  words,  the  component  being  tested  has 
no  "memory"  as  to  previous  stresses  of  a  differ- 
ent type. 

Computation  of  the  system  reliability  is  accomplished 
by  the  following  procedure: 

An  unbiased  failure  rate  estimate  is  computed,  pro- 
viding the  failure  rates  are  anticipated  as  being  small. 
For  ease  of  computation,  a  simplified  unbiased  estimator, 


ath 
ponent,  the  unbiased  estimator,  A,.  ,  is  given  by: 


which  is  slightly  conservative,  is  used.   For  the  lzn   com- 

A 


•*Ibld.   Po3-3. 
**Ibid. 

5ibid. 
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A      fi        2Ni  /  v 

\  =  ir-1 —  * —  (1) 

1    Nj_ 

L        Tn    2N  +  1 

th 
where  N  =  the  sample  size  of  the  i   component  (i.e.  the 

number  of  components  of  type  i  tested). 

T  =  the  test  time  (time  to  failure)  of  the  component 
of  type  i. 

T  .  =  the  planned  test  time  for  component  i. 

and  f  =  the  number  of  RL  failure  times  which  are  less 
i  i 

than  TQi. 

The  system  failure  rate  estimate  is  then  given  by: 


(2) 


where  m  is  the  number  of  components  in  the  system.   The 
component  number  assumes  values  i  =  l,c..,m. 

The  variance  of  the  unbiased  failure  rate  estimate  is 
given  by: 


A. 

m 

A 

\    = 

S 

h 

i=l 

Variance  =  £   U,/S.  ) 
1=1    x   i 


where 


i 

S  =  £     T1 1  =  the  sum  of  all  test  times 
1   j-l   1J 


accumulated  on  the  N  components  of  type  i. 
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These  estimates  of  failure  rates  for  components  In  a 
given  test  condition  can  be  continually  updated  during  the 
development  phase  as  more  data  becomes  available.   Also, 
during  any  higher  level  of  assemblage,  further  testing  can 
also  yield  data  to  modify  the  failure  rate  estimate.   Data 
is  obtained  indicating  the  time  that  each  component  operates 
during  the  system  test. 

The  reliability  estimate  for  the  system  is  then  given 

A 
A        _  \  A 

by  R  =  e   where  k   is  the  appropriate  failure  rate  esti- 
mate  and  only  series  systems  or  subsystems  are  considered, 
given  by  formula  (2). 

The  reliability  estimates  obtained  are  subject  to 
statistical  uncertainty,  therefore  the  important  item  to 
observe  is  the  confidence  interval  about  the  estimate. 
Here,  the  lower  confidence  limit  on  system  reliability 
becomes  the  pertinent  quantity  to  observe.   The  "Guide 
Manual"  bases  the  lower  confidence  limit  on  Normal  theory 
and  it  is  corrected  to  compensate  for  small  values  of  A.. 
The  upper  limit  on  failure  rate,  X    ,  provides  the  corres- 

u  A 

A       _\ 

ponding  lower  confidence  limit  on  reliability  by  R  T  =  e  u. 

S  ,  L» 

The  upper  limit  on  failure  rate  is  given  by: 


^a    2/\     \   .  a  2a    4a2  m 

2\   +  K  c  +   V  4A.K  c  +  K  c 


for  2  f4  >  0 
2  1=1  1 

(3) 

2  m  m 

—  £   (1/S, )  for  £     f .  =  0 


v. 


n  1=1      *•  1=1   i 
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where 


m    .a 


2    U</S.  ) 

C  =  5 


n  =  the  number  of  component-environment-test  condition 
terms  in  the  summation,  and  K  =  the  percentage  point  of  the 
Normal  Distribution.   The  values  of  K  for  a  given  confidence 
level  are  not  completely  appropriate  for  small  values  of  \e 
Consequently,  a  correction  factor,  3,  is  used  to  obtain 
desired  precision.   Beta  values  are  tabulated  for  the  80$ 
confidence  limit  in  reference  (2).   Thus  the  percentage  point 
of  the  Normal  Distribution,  K,  is  modified  by  £K  and  replaced 
in  formula  (3). 

A  more  detailed  description  of  this  reliability  meas- 
urement method  plus  examples  of  the  procedure  are  given  in 
reference  (2). 


m- 


CHAPTER  III 
SIMULATION  PROCEDURE 

Suppose  the  assumption  of  constant  failure  rate  used 
in  the  previously  described  statistical  model  is  relaxedo 
It  is  noitf  desired  to  test  the  accuracy  of  that  model  when 
failures  occur  according  to  a  ^eibull  Failure  Law  and  a  Log 
Normal  Failure  Law.   For  testing  purposes  by  simulation 
techniques,  consider  a  series  system  with  four  components, 
each  component  having  a  well  defined  failure  law. 

The  true  reliability  for  each  component  of  this  series 
system  is  given  by  R  ,  i  =  1,2,3,^,  and  thus  the  true  system 
reliability  is 

R  =   2   R   .  (*0 

s   i=l   i 

th 
Using  the  statistical  model,  an  estimate  of  the  (1-a.) 

a 

percentile  confidence  limit  of  the  system,  R  T,    N,  can  be 

obtained  by  computer  simulation.   This  is  a  random  variable 
as  calculated  from  the  statistical  model.   If 

then,  in  fact,  R  _,  ,  is  an  exact  (1-a)   percentile  lower 

S  ,  J_i  \Q.  ) 

confidence  limit  for  Rs.   This  says  that  Rg  is  always  the 
(1-a)"1  percentile  point  of  the  probability  distribution  of 

A  /\ 

Rs  T(a).   Tlie  distribution  of  R  L(a)  is  constructed  by 
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th 
computer  simulation.   Letting  A-,     be  the  (1-a)   per- 


A 


centlle  point  of  the  distribution  of  R  T ,  * ,  then  the 
absolute  difference,   A-^  _   -  R   ,  is  a  measure  of  the 
accuracy  of  the  statistical  model. 

Two  other  measures  of  the  accuracy  of  this  model  are 
the  estimated  mean,  R  T(a)j  an<^  the  estimated  standard 

deviation,  s£      ,  which  give  some  assurance  that  the 

tts,L(a) 

actual  values  generated  by  the  procedure  are  reasonable. 

The  simulation  procedure  itself  is  now  discussed  for 
the  series  system  of  four  components: 

The  distribution  of  Rg  j-(a\  is  constructed  by  gener- 
ating 500  random  observations  on  R  T(a)  f°r  &  given  system, 
number  of  components  tested,  and  reliability  of  the  com- 
ponents.  For  this  study  a  four  component  series  system  is 
used,  with  each  component  having  reliability  Ri ,  i  =  1,2,3,^. 
The  reliability,  R, ,  is  defined  for  each  component  by 
R.  =  R. (1)  =  P(  T.  >  1) ,  where  T.  is  the  time  to  failure 
random  variable  for  a  component  of  type  i.   This  establishes 
the  parameter  values  for  a  particular  failure  rate  distri- 
bution.  For  example,  if  a  component  is  to  obey  a  Weibull 
Failure  Law, 

R(t)  =  e^^  (5) 

and  R(l)  =  e  A  =  .995,  then,  for  3  =  2,  this  implies  that 
\   =  .0707. 

Similarly,  for  the  Log  Normal  Failure  Law,  if 
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R(t)  =  eZ  =  1  -  §   ln  t   "^  (6) 

and  R(l)  =  1  -  <|(-/V<r  )  =  .995,  then,  for  cr  =  1,  this 
implies  that  yU  =  2.576.   The  failure  rate  functions  deter- 
mined from  the  reliabilities  and  selected  parameters  are 
displayed  in  Appendix  I. 

A  random  number  generator  is  used  to  obtain  a  sample 
of  N.  tested  components  of  type  i  which  fail  according  to 
their  respective  failure  rate  distributions.   This  Monte 
Carlo  method  draws  a  random  number  which  is  Uniformly  Dis- 
tributed (0,1)  and  converts  it  to  a  tfeibull  variate  or  a 

th 
Log  Normal  variate  to  obtain  the  time  to  failure  of  the  i 

component,  T  .   For  the  sA/eibull  Distribution,  e  vai'  =  Y» 

where  Y  is  distributed  as  Uniform  (0,1),  Implies  that 

T  =  j-W^  (7) 

is  a  random  Weibull  variate.   For  the  Log  Normal  Distri- 
bution, if  X  is  a  random  variable  distributed  as  Normal (0,1), 
then  Z  =  cX  +  M.  is  a  random  variable  distributed  as  Normal 
( yA,G-2)  and  thus  T  =  ez  is  a  random  Log  Normal  variate.   In 
this  manner,  N.  failure  times  are  generated  for  each  of  the 
four  components : 


Tll 

T12 

T13  °  ' 

1N]_ 

T21 

T22 

T23  • 

'  "  T2N2 

T3l 

T32 

T33  ' 

'  '  T3N3 

*JH 

T^2 

T^3  • 

'  "  T^% 
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Using  these  failure  times  in  the  statistical  procedure 

A        A 

discussed  in  the  previous  chapter,  values  of  R  and  R_  T(n\ 

S         S  i  J_i  ^OC  ) 

are  obtained  with  that  model.   Each  Ti *    is  examined  for 
failure  before  the  end  of  the  planned  test  time  or  else 
terminated  at  the  end  of  the  planned  test  time.   This  same 
procedure  is  completed  500  times  to  generate  500  random 

A 

observations,  thus  constructing  the  distributions  of  Rs 

A 

and  R     T ,  v. 
s  ,L(ct ; 

The  mean  and  standard  deviation  of  the  distributions 
are  determined  by: 


§._T.t*\  (8) 


A 

Rs,L(cc)   = 

1 
500 

S 

k=l 

R                 1 

s  '    500 

500 

k=l 

A 

ls,L(a)k 


(9) 


\L(a)l5iFl{tS^V!^))2     (10) 


(11) 


The  500  values  of  §g  wa\  are  then  ordered  and  the 

th 
(1-a)   percentile  point  of  the  distribution  is  obtained. 

Thus,  the  value  A-,     is  available  for  comparison  with  R  . 

o-t- 

This  process  is  done  merely  by  counting  down  to  the  101 

A 

value  of  R    ,  .  from  the  largest,  which  represents  the 

S  ,  Xj  \CL  ) 

80th  percentile  point. 
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Several  cases  are  examined  with  various  failure  rate 

distributions  and  with  various  combinations  of  components. 

For  each  case  examined,  the  planned  test  time,  T  .,  varies 

to  produce  an  optimistic  or  pessimistic  reliability  test. 

If  Toi  =  ,5>  the  system  will  yield  an  optimistic  test  since 

the  planned  test  time  will  be  reached  before  most  components 

fail.   The  opposite  will  occur  when  T4  =5.   AT.  chosen 

oi         oi 

such  that  the  average  failure  rate  is  equal  to  the  failure 
rate  at  time  one  appears  to  be  a  good  test  of  system  relia- 
bility for  this  procedure.   That  is,  make  TQ.    such  that  it 
satisfies  the  equation 


z1(t)dt  =  z1(l)  (12) 


where  z(t)  is  the  failure  rate  function;  the  failure  distri- 
bution divided  by  the  reliability  function:  z(t ) =  f (t )/R(t ). 
All  cases  studied  by  simulation  procedures  are  given  in 
Table  I.  All  simulation  results  of  these  cases  are  given  in 
Table  II. 
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TABLE  I 


CASES  STUDIED  BY  SIMULATION  PROCEDURES 


Case 

No. 

Compo- 
nent No. 

Failure 
Law 

Parameters 

Ri 

R 
s 

1 

1 

Wei bull 

3=1.5,  X=.0293 

.995 

.98 

2 

Wei bull 

3=2.0,  X=.0707 

.995 

3 

Log  Normal 

C"=1.0,/l=2.576 

.995 

k 

Log  Normal 

<r  =2.0,/^=5.152 

.995 

2 

1 

Wei bull 

3=1.5,  \=.0293 

.995 

.98 

2 

Wei bull 

3=1.5,  x=.0293 

.995 

3 

Log  Normal 

<r =2.0,^=5.152 

.995 

4 

Log  Normal 

cr=2.o,yA=5„i52 

.995 

3 

1 

Weibull 

3=1.33,  X=.0378 

.987 

.95 

2 

Wei bull 

3=1.33,  \=.0378 

.987 

3 

Log  Normal 

cr=i.5,/^=3.345 

.987 

fc 

Log  Normal 

^=1.5,^=3.3^5 

.987 

** 

l 

Weibull 

3=1.2,  \=.  0^85 

o97^ 

.90 

2 

Weibull 

3=1.2,  ^=.0^85 

.97^ 

3 

Log  Normal 

(1=1.5,^=2.919 

.97^ 

4 

Log  Normal 

(1=1.5,^=2.919 

.97^ 

5 

1 

Weibull 

3=1.2,  \=.  0^+85 

•  97l*> 

.90 

2 

Weibull 

3=1.2,  \=.0^85 

.97^ 

3 

Weibull 

3=1.2,  \=.  0^+85 

.97^ 

4 

Weibull 

3=1.2,  \=.0485 

.97^ 
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TABLE   I    (Continued) 


Case 

No. 

Compo- 
nent No. 

Failure 
Law 

Parameters 

Ri 

Rs 

6 

1 

Log  Normal 

(7=1.5,  >a=2.919 

.974 

.90 

2 

Log  Normal 

<r=l.  5,^=2.919 

.974 

3 

Log  Normal 

a-=1.5jyu=2.919 

.974 

4 

Log  Normal 

cr=1.5,A=2.9l9 

.974 

7 

1 

deibull 

3=1.1,  \=o0^83 

.965 

.867 

2 

tfeibull 

3=1.1,  x=.0483 

.965 

3 

rfeibull 

(3=1.1,  \=.0^83 

.965 

4 

v/eibull 

3=1.1,  \=.0483 

.965 

8 

1 

Log  Normal 

C=2.0,  ^=4.0 

.977 

.912 

2 

Log  Normal 

(T=2.  0,^=^.0 

.977 

3 

Log  Normal 

(T=2.0,/a.=^.0 

.977 

** 

Log  Normal 

<T=2.0,  >u=i+.0 

.977 
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CHAPTER  IV 
RESULTS  AND  CONCLUSIONS 

In  the  preceding  chapters  the  statistical  procedure 
from  the  Guide  Manual  For  Reliability  Measurement  program 
has  been  briefly  explained  to  obtain  a  100(1  -  a)%   lower  con- 
fidence limit  on  system  reliability.   The  simulation  of  this 
method  on  a  computer  has  been  explained  and  the  technique 
for  measuring  the  accuracy  of  the  statistical  procedure 
stated.   Now,  from  the  results  of  the  simulation,  some  con- 
clusions may  be  drawn. 

The  accuracy  of  the  statistical  procedure  has  been 
examined  for  eight  cases  which  represent  different  combi- 
nations of  components  for  various  sets  of  parameter  values. 
The  results  for  these  cases  are  given  in  Table  II.   The 
mean  and  standard  deviation  of  the  500  values  of  R  T(a\ 
and  Rg  are  given,  along  with  the  (l-a)^  percentile  point, 
A-j_  _  a.   In  all  cases  studied,  a  was  taken  to  be  .20. 

The  accuracy  of  the  procedure  is  measured  by  the 
absolute  difference  between  the  80   percentile  of  the  dis- 

A 

tribution  of  R        ,    20P   whiclrl  is   A  80'    anci  the  actual   system 


reliability,  R_  (i.e.  accuracy  = 


A.80  "  Rs 


Also  stated  in  the  results,  is  the  quantity  TT,  which 
is  the  amount  of  testing  relative  to  the  component  un- 
reliabilities.  In  other  words,  this  is  a  measure  of  the 
amount  of  testing  required  to  achieve  a  desired  accuracy 
for  a  given  system  with  component  failure  rates,  z,(l), 
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and  planned  test  times,  T  .  .   Then  TT  is  given  by 


m 


TT  =   £   N  z  (1)T    .  (13) 

i=l  X    x     01 


As  an  example  of  the  accuracy  study,  consider  Case  1, 
where  the  system  consists  of  four  components  in  series  dis- 
playing the  following  failure  laws: 


Component  1 
Component  2 
Component  3 
Component  k 


^Teibull,  pi  =  1.5,  ^=.0293 
ieibull,  (3  =  2.0,  k=.C?07 
Log  Normal,  G"  =  1.0  ,  yu.  =  2.576 
Log  Normal,  <r  =  2.0,  yiA=  5.152 


The  component  reliabilities,  R,  =  .995,  for  all  four 

components,  so  the  system  reliability,  R  =  . 9801.   The 

s 

measure  of  accuracy  for  the  optimistic  planned  test  time, 

I 

I  .  =  .5»  i8  A  oA  -  L   =  .008  for  500  components  of  each 
oi  1  .00    s 

type  tested.   The  pessimistic  planned  test  time,  TQi  =  5.0, 
yields  an  accuracy  of  .055  for  100  components  of  each  type 
tested.   The  ad  hoc  planned  test  times,  T  ,  =  2.25,  To2  =  2.0, 
TQo  =  1.95,  and  T  .,  =  3.25,  as  determined  from  formula  (12), 
give  an  accuracy  of  .030  for  !L  =  100.   The  accuracy  for 
the  ad  hoc  planned  test  times  is  only  fair.   This  planned 
test  time  represents  the  best  measure  of  accuracy  of  the 
three  selected,, 

It  can  be  seen  from  Table  II  and  Appendix  I  that  the 
accuracy  of  the  procedure  is  a  function  of  how  closely  the 
tfeibull  and  Log  Normal  failure  rate  functions  approximate 
a  constant  failure  rate  function.   The  greater  the  t/Jeibull 
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and  Log  Normal  failure  rate  functions  deviate  from  the 
constant  failure  rate,  the  more  inaccurate  the  procedure. 
Also,  the  accuracy  decreases  as  the  planned  test  time  is 
increased.   The  Weibull  failure  rate  function  increases  in 
time  and  the  Log  Normal  failure  rate  function  increases 
then  decreases  in  time,  thus  deviating  from  the  constant 
failure  rate  for  large  planned  test  times. 

It  can  thus  be  concluded  from  the  study  that  the 
statistical  procedure  can  be  useful  for  reliability  approx- 
imations of  a  system,  if  the  components  that  fail  according 
to  these  non-constant  rates  do  not  deviate  too  significantly 
from  a  constant  failure  rate. 
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APPENDIX  I 


GRAPHS  OF  COMPONENT  FAILURE  RATE  FUNCTIONS 
FOR  EACH  SYSTEM  TESTED 
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FIGURE  1 
Component  Failure  Rate  Functions,  z(t),  for  Case  1. 
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FIGURE  2 
Component  Failure  Rate  Functions,  z(t),  for  Case  2. 
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APPENDIX  I   (Continued) 
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FIGURE  3 
Component  Failure  Rate  Functions ,  z(t ) ,  for  Case  J, 
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APPENDIX  I   (Continued) 
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FIGURE  k 
Component  Failure  Rate  Functions,  z(t),  for  Case  4. 
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APPENDIX  I   (Continued) 
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FIGURE  5 
Component  Failure  Rate  Functions,  z(t),  for  Case  5. 
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APPENDIX  I   (Continued) 
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FIGURE  6 
Component  Failure  Rate  Functions,  z(t),  for  Case  6- 
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APPENDIX  I   (Continued) 
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FIGURE  7 
Component  Failure  Rate  Functions,  z(t),  for  Case  7. 
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FIGURE  8 
Component  Failure  Rate  Functions,  z(t),  for  Case  8, 
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APPENDIX  II 


COMPUTE*  PROGRAM 


DIMENSION  PS(5OC),RSL(50C) ,S(4) ,NF(4I, 
*N( A) ,HLAM( A) tRTA(30) ,Tn(4,1) 

C      RFAO  IN  SAMPLE  SI7E,  N 
READ  1 ,M 

1  FORMAT  (418) 

C      READ  IN  BFTA  CORRECTION  FACTORS 
REAO  2,  RTA 

2  FORMAT  (5FR.3) 

C      READ  IN  PLANNED  TEST  TIMES 

RCAD  3,  UTOU,J>,I»l,4l, J*l, 31 

3  FORMAT  (4F10.0) 

C      RFAD  IN  FAILURE  DISTRIBUTION  PAR4MFTFPS 

READ  12,  ALFA1,BETA1, ALFA2,BFTA2, 
♦SIGMA! ,XMFAN1 ♦ SIGMA2, XMEAN2 
1 2  FORMAT  (8F6.0) 

PRINT  4 
A  PORMAT  (1H1 ) 
C      RUN  SIMULATION  FHR  THREE  OIFFEPENT  PLANNED 
C      TEST  TT^FS 

DO  <399  TPUN  =  1,3 
C      URN  IS  RANDCM  NUMBER  GENERATOR 

Y  =  1JRMO) 

RTOT=C. 

RLTOT=D. 
C  GFNERATE    500    VAL'tFS    nc    LOW=R    CONFIDENCE    LIMIT 

DO    SO    K  =  1  ,*0n 

HLAMD=^. 

FTOT=0. 

TS  =  0. 

TSUM=n. 
C       SIMULATE  FAILURE  OF  FACH  OF  THE  FOUR  COMPONENTS 

DO  30  1=1,4 

S(I )=^. 

NF(I)=0 

*=N(I) 
r  GENERATE    N    FAILURE    TIMES    Fnp    FACH    COMPONFMT 

DO    2^    J=1,m 

r?n    TO    (5  ,*,7,R)  ,  T 
c    Y=URN(1) 

BTNV1-1./BETA1 

TP=(  (-ALOG(  Y)  )**RINVU/ALFAl 

GO  TO  1C 

6  Y=URM(T ) 
«MNV2=1./BETA2 

TP=(  (-ALOGI  V)  )**BINV?)/ALFA2 
GO  TO  10 

7  TSl=^. 

DO  9  IJ=1,12 
C  TS1=TS1*URN(1 ) 
X!=T<:!-6. 

71=SIGMA1*X1*XMEAN1 
TP=EXP(7\) 
GO  TO  10 
P  TS2=o. 

DO  11  IJ=1,12 
11  TS2=TS2+URN(1 > 
X2=TS2-6. 

Z2=SIGMA2*X2^XMEAN? 
TP=EXP(Z2) 
10  CONTINUE 
C      COUNT  NUMBER  OF  FAILURES  BEFOPF  PLANNED 
C       TEST  TIMP  FXPIRFS 

IF(TP.LT.TO(I,TRUN) >  GOTO  15 
T=TO(T , IRUN) 
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IB 
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32 
3? 
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41 
5^ 


6C 


61 


65 
70 

pr 


or 

qcc 


an  td  18 

T  =  TP 

NF( I)=NF(I  ) 

S(  T  )  =  S(  I  UT 

CONTTMUE 

CALCULATE  L 

RFLIARILITY 

XN=N( I ) 

E  =  NF( T ) 

HLAMH)  =  (F/ 

ftot=ftot«-f 
hlamd=hlamo 
tsum  =tsum«- 

TS  =  TS«-1, /S( 

CONTINUE 
RS(K)=EXP(- 
RTOT=RTOT*R 
CALCULATE  (I 
LOWER  CONE  I 
IPT0T=FT0T4- 
TF(  IETOT-3' 
8ETA=BTA(IF 

co  to  33 

RFTA=1. 
X«=(.R42>*R 
IMETOT)  3* 
ULAMD=( XK** 
GO  TO  41 
CHAT=TSUM/H 
ULAMD=(2, *H 
*( XK**?)*CHA 
RSL(K)=EXP( 
R|.TOT  =  RLTGT 

CONTINUE 

CALCULATE  M 
PR AR=RTOT/* 
SU*=0. 

no  60  K=l,5 

OIFF=(RS(K) 
SI|H  =  SUM-fDlF 

CONTINUE 
SR=S0OT(SUM 
RLRAP=RLTHT 
Sijm=o. 

00  61  *  =  1,5 
r>IPF  =  (RSL(K 
SUM  =  SlJM*niF 
CONTINUE 
SRL=SQRT(SU 
ORDER  500  V 
ON  RFLTABIL 
11  =  1 

no  so  k  =  i,i 

B=RSL(1) 

00  70  1=1,500 

)  65,70,70 


♦  1 


AMnA  HAT  AND  SIMULATED  SYSTEM 


S(I )»*(  2.  +  XN/(2.*XN-H.  )  ) 

♦HLAM(I » 

(HLAM(I )/5(I|| 
II 

HL Awni 

S<K) 

PPER  LAMDA  HAT  AND 

DFNCF  LIMIT  RFLIARILITY 

1- 

)  31,32,32 

TOT) 


ETA 

2)*TS/4. 

LAMO 

LAM0-MXK**?)*CHAT  +  SQRT(4.*HL  AMD* 

T+(XK**4)* (CHAT**2) ) )/? , 
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APPENDIX  III 

c 

WEIBULL  DISTRIBUTION 

The  v/eibull  Probability  Distribution  is  defined  by: 

Q 

f(t)  =  pxM"V  ut)  o<t<  oo 

where  \   is  the  "scale"  parameter,  and  3  is  the  "shape" 
parameter. 

The  Cummulative  Distribution  Function  (CDF)  is  then 


f(x)dx  =  1  -  e  ~  Ut) 


and  the  failure  rate  function  is 

z(t)  =  f(t)/R(t)  =  \^^  _1 

where  R(t)  =  1  -  F(t)  =  the  reliability  function.   The 
Jeibull  Probability  Distribution  reduces  to  the  Exponential 
Probability  Distribution  when  3=1.   Increasing  failure 
rates  are  obtained  only  when  3  >  1. 


"Daniel  R.  Cox,  Renewal  Theory  (London;  Methuen  and 
Co.  Ltd. ,  1962)  p.  21. 
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APPENDIX  IV 


LOG  NORMAL  DISTRIBUTION 


7 


The  Log  Normal  Probability  Distribution  is  given  by 


f(t)  = 


1_ 

2 


In  t  -jo. 


VZTT^t 

=  e* 


0  <  t  <  oo 


where  X  is  distributed  as  Normal  (/jl9g~  ).  For  computational 
purposes,  it  is  beneficial  to  derive  the  Cummulative  Distri- 
bution Function  as 

F(t)  =  p(  t<T  ) 
=  p(  t  < ex  ) 
=  p(ln  t  <  X  ) 


=  P 


In    t    -yW 

<  z 

(T 

In   t    -.M 

G- 

then 


-ft 


f(t)  =  F'(t)  =  -±-   (|) 


In  t  -jm 


The  failure  rate  function  is  given  by 


z(t)  =  f(t)/R(t)  = 


crt   L 


In  t  -  M 


(T 


-t 


In  t  -/a 


7  Paul  L.  Meyer,  Introductory  probability  and  Sta- 
tistical Applications  (Reading,  Mass:  Addison- Lesley,  19&5) 
p.  187. 
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